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We study the explosion mechanism of collapse-driven supernovae by numerical sim- 
ulations with a new nuclear EOS based on unstable nuclei. We report new results of 
simulations of general relativistic hydrodynamics together with the Boltzmann neutrino- 
transport in spherical symmetry. We adopt the new data set of relativistic EOS and the 
conventional set of EOS (Lattimer-Swesty EOS) to examine the influence on dynamics 
of core-collapse, bounce and shock propagation. We follow the behavior of stalled shock 
more than 500 ms after the bounce and compare the evolutions of supernova core. 



1. Introduction 

Understanding the explosion mechanism of core-collapse supernovae is a challenging 
problem, that requires extensive researches in nuclear physics and astrophysics. In order 
to reach the final answer, it is necessary to investigate the core-collapse supernovae by 
implementing hydrodynamics and neutrino-transfer together with reliable nuclear equa- 
tion of states and neutrino-related reactions. In this regard, recent numerical simulations 
of neutrino-transfer hydrodynamics |l|2|3|4j have cast a light to the importance of nu- 
clear physics as well as neutrino-transfer, though they show no explosion at the moment. 
Meanwhile, advances in physics of unstable nuclei have given chances to provide us with 
nuclear physics in supernovae than ever before, therefore, those new nuclear data should 
be examined in modern supernova simulations. In this paper, we focus on the influence 
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of the new nuclear equation of state (EOS) in the neutrino-transfer hydrodynamics. We 
follow the core-collapse, bounce and shock propagation by adopting the new equation of 
state, which is based on the data of unstable nuclei, and the conventional one which has 
been used almost uniquely in recent simulations. We compare the behavior of shock and 
the thermal evolution of supernova core by performing numerical simulations for a long 
period of ~1 sec after the core bounce. 

2. A new nuclear EOS table 

Recently, a new complete set of EOS for supernova simulations (Shen's EOS) has be- 
come available jolGj . The relativistic mean field (RMF) theory with a local density approx- 
imation has been applied to the derivation of the supernova EOS table. The RMF theory 
has been successful to reproduce the saturation properties, masses and radii of nuclei, and 
proton-nucleus scattering data [7j. The effective interaction used in the RMF theory is 
checked by the recent experimental data of unstable nuclei in neutron-rich environment 
close to astrophysical situations [EE]- 

We stress that the RMF theory jHj is based on the relativistic Briickner-Hartree-Fock 
(RBHF) theory |lUj . The RBHF theory, which is a microscopic and relativistic many 
body theory, has been shown to be successful to reproduce the saturation of nuclear matter 
starting from the nucleon-nucleon interactions determined by scattering experiments. This 
is in good contrast with non-relativistic many body frameworks which can account for 
the saturation only with the introduction of extra three-body interactions. 

The RMF framework with the parameter set TM1, which was determined as the best 
one to reproduce the properties of finite nuclei including n-rich ones jS], provides the 
uniform nuclear matter with the incompressibility of 281 MeV and the symmetry energy 
of 36.9 MeV. The maximum neutron star mass calculated for the cold neutron star matter 
in the RMF with TM1 is 2.2 M [TTJ . The table of EOS covers the wide range of density, 
electron fraction and temperature, which is necessary for supernova simulations. The 
relativistic EOS table has been applied to numerical simulations of r-process in neutrino- 
driven winds ^2] and prompt supernova explosions [13*] . and other simulations |14|15|16j . 

For comparison, we adopt also the EOS by Lattimer and Swesty J7j. The EOS is 
based on the compressible liquid drop model for nuclei together with dripped nucleons. 
The bulk energy of nuclear matter is expressed in terms of density, proton fraction and 
temperature with nuclear parameters. The values of nuclear parameters are chosen to 
be the ones suggested from nuclear mass formulae and other theoretical studies with the 
Skyrme interaction. Among various parameters, the symmetry energy is set to be 29.3 
MeV, which is smaller than the value in the relativistic EOS. As for the incompressibility, 
we use the EOS with 180 MeV, which has been often used for recent supernova simulations. 

3. Simulations of core-collapse supernovae 

We have developed a new numerical code of neutrino-transfer hydrodynamics |18|19|2U] 
for supernova simulations. The code solves hydrodynamics and neutrino-transfer at once 
in general relativity under the spherical symmetry. The Boltzmann equation for neutrinos 
is solved by a finite difference scheme (Sat method) together with lagrangian hydrodynam- 
ics in the implicit manner. The implicit method enables us to have a longer time step than 
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the explicit method, therefore, our code is advantageous to follow long term behaviors 
after the core bounce. We adopt 127 spatial zones and discretize the neutrino distribution 
function with 6 angle zones and 14 energy zones for u e , u e , v^/ T and p^/ti respectively. 
The weak interaction rates regarding neutrinos are based on the standard rates by Bruenn 
[2~Tj . In addition to the Bruenn's standard neutrino processes, the plasmon process and 
the nucleon-nucleon bremsstralung processes are included [20] • The collision term of the 
Boltzmann equation is explicitly calculated as functions of neutrino angles and energies. 
As an initial model, we adopt the profile of iron core of a 15M progenitor from Woosley 
and Weaver |22,. We perform the numerical simulations with Shen's EOS (denoted by 
SH) and Lattimer-Swesty EOS (LS) for comparisons. 
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Figure 1. Radial trajectories of mass ele- 
ments of the core of 15M star clS db function 
of time after the bounce in SH. 
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Figure 2. Profiles of entropy, temperature 
and lepton fraction at 600 ms after the 
bounce in SH (thick) and LS (thin). 



It is remarkable that the explosion does not occur in the model SH, i.e. the case with 
the new EOS table, and the shock wave stalls in a similar manner as the model LS. We 
have found, however, that there are quantitative differences in core-collapse and bounce 
due to the differences between nuclear EOSs. The peak central density at the bounce 
in the case of SH is 3.3 xlO 14 g/cm 3 , which is lower than 4.2 xlO 14 g/cm 3 in the case of 
LS. It is to be noted that, in microscopic many body calculations, the relativistic EOS 
(including Shen's EOS) is generally stiffer |TU] than the non-relativistic EOS, on which 
Lattimer-Swesty EOS is based. The stiff EOS is not advantageous in a simple argument 
of initial shock energy, however, the symmetry energy of EOS also plays an important 
role ■ Having a larger symmetry energy, free-proton fractions during collapse in SH 
are smaller than in LS. Smaller free-proton fractions lead to smaller electron capture rates 
when electron captures on nuclei are suppressed. Because of this difference, the trapped 
lepton fraction at the bounce in SH is 0.35 at center, which is slightly larger than 0.34 in 
LS. In the current simulations, this difference of lepton fraction does not change the size 
of bounce core significantly. 
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As a result, the shock in SH does not reach at a significantly larger radius than in LS 
(Fig. P). The shock stalls below 200 km and starts receding in two cases. The central core 
becomes a proto-neutron star having a radius of several tens km and a steady accretion 
shock is formed. The difference of EOS becomes apparent having a more compact star for 
LS with the central density of 6.0 xlO 14 g/cm 3 (3.9 xlO 14 g/cm 3 for SH) at 600 ms after 
the bounce. The peak structure of temperature at around 10 km is formed and the peak 
temperature reaches at ~40 and ~50 MeV for SH and LS, respectively, due to the gradual 
compression of core having accretion (Fig. |2J). It is noticeable that negative gradients in 
entropy and lepton fraction appear by this stage, suggesting the importance of convection. 
These differences of thermal structure may give influences on shock dynamics, supernova 
neutrinos and proto-neutron star cooling. Further details of the full numerical simulations 
will be published elsewhere |2T)] . 
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